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Abstract 

The problem of radio wave reflection from an opticaUy thick plane monotonous 
layer of magnetized plasma is considered at present work. The plasma electron 
density irregularities are described by spatial spectrum of an arbitrary form. 
The small-angle scattering approximation in the invariant ray coordinates is 
suggested for analytical investigation of the radiation transfer equation. The 
approximated solution describing spatial-and-angular distribution of radia- 
tion reflected from a plasma layer is obtained. The obtained solution can be 
applied, for example, to the ionospheric radio wave propagation. 

1 Introduction 

Basic goal of the present work consists in derivation of the transfer equation solution 
describing spatial-and-angular distribution P (p, lu) of radio radiation reflected from 
a plane stratified layer of magnetized plasma with random irregularities. 

The radiation transfer equation (RTE) in a randomly irregular magnetized 
plasma was obtained in the work [1] under rather general initial assumptions. In 
particular, the medium average properties were assumed smoothly varying both in 
space and in time. In the work [2] the radiation energy balance (REB) equation 
describing radiation transfer in a plane stratified layer of stationary plasma with 
random irregularities has been deduced. The invariant ray coordinates, allowing 
one to take into account by a natural way refraction of waves and to represent the 
equation in the most simple form, were used there. In the work [3] it was shown that 
the equation REB is a particular case of the radiation transfer equation obtained 
in [1] and can be deduced from the latter by means of transition to tlie invariant 
ray coordinates. Equation REB, thus, allows one to investigate influence of multiple 
scattering in a plane stratified plasma layer on the characteristics of radiation. In 
particular, it enables one to determine the spatial-and-angular distribution of ra- 
diation leaving the layer if the source directivity diagram and irregularity spatial 
spectrum are known. A few effects which require of wave amplitudes coherent sum- 
mation for their description (for example, phenomenon of enhanced backscattering) 
are excluded from consideration. However, the multiple scattering effects are much 
stronger, as a rule. This is particularly true for the ionospheric radio propagation. 

The numerical methods of the transfer equation solving developed in the theory 
of neutron transfer and in the atmospheric optics appear useless for the equation 
REB analysis. They are adapted, basically, to the solution of one-dimensional prob- 
lems with isotropic scattering and plane incident wave. In a case of magnetized 
plasma the presence of regular refraction, aspect-sensitive character of scattering 
on anisometric irregularities and high dimension of the equation REB (it contains 
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two angular and two spatial coordinates as independent variables) complicate con- 
struction of the effective numerical algorithm for its solving. In this situation it 
is expedient to solve the equation REB in two stages. The first stage consists of 
obtaining of the approximated analytical solution allowing one to carry out the 
qualitative analysis of its properties and to reveal of its peculiarities. At the second 
stage the numerical estimation methods can be applied to the obtained analytical 
solution, or methods of the numerical solving of the initial equation taking into 
account the information obtained at the first stage can be designed. Therefore the 
problem of obtaining of the equation REB approximated analytical solutions is of 
interest. 

We begin the present paper from a detailed exposition of the invariant ray 
coordinates concept. Then possibility to use of the small- angle scattering in the 
invariant coordinates approximation is discussed. Two modifications of the REB 
equation solution are obtained. The analysis of the obtained solutions concludes 
the paper. 

2 Invariant ray coordinates and the radiation en- 
ergy balance equation 

It is convenient to display graphically the electromagnetic wave propagation in 
a plane-stratified plasma layer with the aid of the Poeverlein construction [4,5]. 
We shall briefly describe it. Let the Cartesian system of coordinates has axis z 
perpendicular and the plane xOj/ parallel to the plasma layer. We shall name such 
coordinate system "vertical" . It is assumed that the vector of the external magnetic 
field H is situated in the plane z^y . Module of radius-vector of any point inside 
of the unit sphere with centrum in the coordinate origin corresponds to the value 
of refractive index ni(u,Q!), where i = 1 relates to the extraordinary wave, i = 
relates to the ordinary one, v = w^/w^ , is the plasma frequency, u)"^ is the 
frequency of a wave, a is the angle between radius-vector and magnetic field H . 
The refractive index surface corresponding to a fixed value of v and to all possible 
directions of the radius- vector represents a rotation body about an axis parallel to 
vector H (see fig. 1). 

Convenience of the described construction (in fact, this is an example of coordi- 
nate system in space of wave vectors k) is become evident when drawing the wave 
trajectory: it is represented by a straight line, parallel to the axis z . This is a 
consequence of the generalized Snell law, which also requires of equality of the fall 
angle and exit angle onto/from a layer {9 ), and constantness of the wave vector 
azimuth angle ( ip). Note, that the crossing point of a wave trajectory with the 
refractive index surface under given value of v determines current direction of the 
wave vector in a layer (it is anti-parallel to a radius-vector) and current direction 
of the group speed vector (it coincides with the normal to the refractive index sur- 
face). The projection of a wave trajectory onto the plane xOy is a point which 
radius- vector has module ainO and its angle with relation to axis x equals to ip . 
Thus, the coordinates define completely the whole ray trajectory shape in a plane 
layer and outside of it and are, in this sense, invariant on this trajectory. 

Radiation of an arbitrary point source of electromagnetic waves within the solid 
angle 6 + d9; ^ ip + dip corresponds to the energy flux in the fc-space inside 
of a cylindrical ray tube parallel to axis z with cross section sm9d{sin9)dip = 
sin 9 cos 9d9dp . In case of regular (without random irregularities) plasma layer 
this energy flux is conserved and completely determined by the source directivity 
diagram: 
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p{z-e,^,p) = p^{e,^,p) , (1) 

where P is energy flux density in the direction determined by angles 9, ip through 
the point p on some base plane situated outside of the layer parallel to it (in the 
ionosphere case it is convenient to choose the Earth' surface as the base plane), z is 
distance from the base plane (height in the ionosphere case). We shall assume in the 
present paper that function z{v) is monotonous in the region of wave propagation 
and reflection. If random irregularities are absent and source of radiation is point, 
variable p in (|l|) is superfluous, as the matter of fact, since unequivocal relation 
between it and angles of arrival of a ray 0, exists. When scattering is present 
the radiation energy redistributes over angular variables 0, 95 and in space what is 
described by variable p . The value of P satisfies in this case to the equation of 
radiation energy balance [2,3]: 
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^^P{z,e,ip,p) = J{-P{z;9,ip,p)sm9coseC-^{z;9,(p) 
■a [ao {9, ip) , /3o (9, ^) ; a (9', if') , /3 [9', ^')] sma{9' , 
+P 



d(a,l3) 



z;e',tp',p~ ^z;9',tp';9,ip) smO' cose'C~^{z;9' ,ip')- 



(2) 

a [ao {9', ^'),(3o {9', v'); a {9, v), (3 {9, f)] sin a{9,v) }d9'd^' 



C{z]9,ip) is cosine of a ray trajectory inclination angle corresponding to the in- 
variant angles 9 and ip;\d{a, (3) I d{9,Lp\ is Jacobean of transition from angular 
coordinates 9, ip to the wave vector polar and azimuth angles a and [3 in the 
"magnetic" coordinate system (which axis O2 is parallel to the magnetic field); 

[ao {9, ^) , /?o {0, f) ; a {9', ip') , /3 {9' .p>')] = a \9, ^; 9', ip'\ 

is scattering differential cross section describing intensity of the scattered wave 
with wave vector coordinates a, /? in magnetic coordinate system (corresponding 
invariant coordinates are 9' and ip') which arises at interaction of the wave with 
wave vector coordinates q;o,/3o (invariant coordinates 9 and p> ) with irregularities. 
Vector function <i>(z; 6*', p' \ 9, p) represents the displacement of the point of arrival 
onto the base plane of a ray which has angular coordinates 9' and p' after scattering 
at level z with relation to the point of arrival of an incident ray with angular 
coordinates 9,p . It is essential that in a plane-stratified medium the function $ 
is determined only by smoothed layer structure v{z) and does not depend on the 
scattering point horizontal coordinate and also on coordinate p of the incident 
and scattered rays. Note also that ratio $(z; 9, p; 9' , p') = —$(2;; 9' , p'; 9, p) takes 
place. 

It is possible to check up that equation (2) satisfies to the energy conservation 
law: when integrating over all possible for level z values of 9, p and all p its right 
side turns into zero. It is natural since in absence of true absorption the energy 
inside the plasma layer does not collected. 

Analyzing expression for the scattering differential cross section in a magnetized 
plasma (see, for example, [6]), it is easy to be convinced that the following symmetry 
ratio takes place: 

cr [9, p; 9', p'\ r? cosiJ^ = a \9' , p'; 9, p] n'^ cosiJg (3) 

where i?g is angle between the wave vector and group speed vector, n is refractive 
index. Using (Q) the equation (^) can be presented as follows: 

f^P{z,p,9,p)^ jQ(z;9,p;9',p') 

Piz,p- <Hz;9',p';9,p>),9',p') ~ Piz,p,9,p)} d9'dp' ^ > 
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where Q{z;6,ip;9f,(ff) — (j{d,Lp;9' ,(p')C ^ {z , 9 , (p) sin 6' \dil'f^ / dil' \ , and quantity 

Q {z;9,(p;9' ,(p') = Q{z;9,ip;9',(p')sm9cos9 is symmetric with relation to rear- 
rangement of shaded and not shaded variables. The equation REB in the form 
has the most compact and perfect appearance. It is clear from physical reasons that 
(||) has to have the unique solution for given initial distribution Po{9, ip, p). The ob- 
tained equation can be directly used for numerical calculation of the signal strength 
spatial distribution in presence of scattering. However, as it was noted at intro- 
duction already, this approach leads to essential difficulties. Subsequent sections 
describe the method of construction of the energy balance equation approximated 
analytical solution. 



3 Small-angle scattering approximation in the in- 
variant ray coordinates 

Let us consider the auxiliary equation of the following kind, which differs from (4) 
only by absence of the dash over variable lo marked by arrow: 



■^P{z,p,uj)^ J Q{z-iu-u;')I^P{zJ+ ${z;iu;iu'),u)-P{z,p,iu)^du;' (5) 

where designation uj = {9, , duo = d9dip has been used for the sake of compact- 
ness. Equation (^) can be easily solved analytically by means of Fourier transfor- 
mation over variable p. The solution has the following form: 

P{z,q,uj) = PQ{q,uj)S{z,0;q,uj), (6) 

where PQ{q,u}) is the Fourier image of the radiation energy flux density passing the 
layer in absence of scattering and the value of S is defined by the expression 



S{z2,zi,q,uj) = exp 



dz' 



du)'Q{z'; uj; u') exp {iq^{z'\ oj; oj'), uj) 



(7) 

One should note that integration over z in this and subsequent formulae, in fact, 
corresponds to integration along the ray trajectory with parameters 0, Lp. The area 
of integration over lo' includes rays which reflection level hr{uj') > z . 
Let us transform now equation (||) by the following way: 



£P{z, p,uj) = J duo'Q{z- UJ- uj') |p(z, p + w; uj'), uj) - P{z, p, w)} 
+ ^ duj'Q{z;uj\uj')^P{z,p+ ^[z-uj;uj'),uj') - P(z, p + ${z]uj;uj'),u})'^ 

Its solution will be looked for in the form 



(8) 



P{z,p,uj)^P {z,p,uj)+X{z,p,uj) (9) 

Thus, auxiliary equation (||) allows to present the solution of the equation (|4|) in 
the form (|). This is an exact representation while some approximated expressions 

for quantities P and X are not used. 

By substituting of (^) into the equation one can obtain the following equation 
for the unknown function X : 
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±X{z, p,Lu)^J duj'Q{z- UJ- uj'){p [z, p+ $(z; lo; lo'), u')- 

- P (z, p + $(z; UJ- uj'), uj)] + / <hj'Q{z-, UJ-, uj'){X{z, p+ ${z; w; uj'), u')- (10) 
-X{z,p,uj)} 

We shall assume now that the most probable distinction of angles uj' and uj is 
small. The heuristic basis for this assumption is given by analysis of the Poeverlein 
construction (fig. 1). It is easy to be convinced examining the Poeverlein construc- 
tion that scattering near the reflection level even for large angles in the wave vectors 
space entails small changes of the invariant angles 9, (p. This is especially true for 
irregularities strongly stretched along the magnetic field (in this case the edges of 
scattered waves wave vectors form circles shown in fig. 1 as patterns A and B). 
One should note also that the changes of invariant angles 6, ip are certainly small 
if scattering with small change of a wave vector direction takes place. This situa- 
tion is typical for irregularity spectra, in which irregularities with scales more than 
sounding wave length dominate. Thus, the small-angle scattering approximation 
in the invariant coordinates has wider applicability area than common small-angle 
scattering approximation. 



Scattering with small changes of 6, ip entails small value of $ 
directly both from sense of this quantity and from the fact what 



That follows 



$(z,u;,w) 







. Let us make use of that to carry out expansion of quantity X at the right side 
of the equation (10) into the Taylor series with small quantities uj' ^ uj and $ . 

Note that making similar expansion of function P at the initial equation (^) would 
be incorrect since function P may not to have property of continuity. For example, 
in case of a point source, Pq is a combination of S -functions. As it will be shown 
later, the function X is expressed from Pq by means of repeated integration and, 
hence, differentiability condition fulfills much easier for it. 

Leaving after expansion only small quantities of the first order, we obtain the 
following equation in partial derivatives: 

d d d 

—X{z,p,uj) - A^{z,uj)—X{z,p,uj) + Af;{z,uj)—X{z,p,uj) = /( z,p,uj), (11) 
where 

Ai^{z,uj) = J duj'Q{z-,uj-,uj'){uj' - uj)-, 
Af;(z,uj)^ / duj'Q{z-,uj-,uj')${z,uj,uj')-, 



/( z, p,uj)^ J duj'Q{z-, UJ-, uj') 

{p{z,p+ ^z;uj;uj'),uj')- P{z,p+ S{z-,uj-,uj'),uj)^ 
Here is the characteristic system for the equation (pj): 



dX J,, ^ ^ dp duj . , . 

-J- = J[ Z:P,^)i T = ^py^^^h -1- = -A^{z,uj), 

az dz dz 

and initial conditions for it at z = 



X ^0-, p = p; = Wo. 
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It is necessary to emphasize the distinction between quantity (! , which is a function 
of z , and invariant variable p. 

Solving the characteristic system we obtain: 

w = a;(z,a;o), ^ = P~ dz'ylp'(z', w(z', wq)) ■ 

J z 

where zq is z -coordinate of the base plane. It follows that 



X{zq,p,lu) 



dz'fl z',p- 



dz"Ap~[z",iu{z",u;o)],io{z',u;o) 



(12) 



Generally, expression (|T^) gives the exact solution of the equation ([TTI). However, 
since we are already within the framework of the invariant coordinate small-angle 
scattering approximation which assumes small value of Ai^(z,lu) , it is possible to 
simplify the problem a little. Assuming = and omitting index at invari- 
ant coordinates to , we are coming to the following approximate representation for 
function X : 



X{zo,p,Lj)^ Jo"dz'{p z,p+ S{z';co;uj') + D{zo,z',Lu),Lu' 



P 



z',p+ <^{z;iLj;iLj') + D{zo, z' ,uj),UJ } 



(13) 



where D{z2, zi,ijj) — J^^ J duj'Q{z;uj;uj')^{z,uj,u!') . 

Thus, in the invariant coordinate small-angle scattering approximation the so- 
lution of the equation REB (||) is represented as a sum of two terms (see the 
first of which is 



P{z,p,U}) = ^ Jd^qPo{q,Ujy 



(14) 



■exp|iqp + /gMz'/ dw'Q(z'; cj; w') exp ^iq<i>(z'; w'), w)^ -1 | 

where J d^QPoi'lT'^) ^^Pi^QP) =^ Po{P:^): a-nd the second one is given by ex- 

pression (p^. 

The solution can be presented in the most simple form if one uses again the 
smallness of quantity $ and expands the second exponent in the formula ( |l^ 
into a series. Leaving after expansion only small quantities of the first order, one 
can obtain: 



P{z,p,uj) = Po p + D{zo,0,uj),uj + J™ dz' J duj'Q{z';uj;uj') 
■{P z,p+ ^{z';uj;uj') + D{zo,z',uj) + D{z',0,uj'),uj' 



(15) 



P 



z',p+ ^{z;uj:uj') + D{zo,Q,uj),uj } 



The last operation is the more precise the faster value of Po{q,uj) decreases under 
\q\ — > oo. The solution of the radiation energy balance equation obtained in 
the present section in the form (^, (p^, (p^), or in the form (p^), expresses the 
spatial-and-angular distribution of radiation intensity passing layer of plasma with 
scattering through the spatial-and-angular distribution of the incident radiation, 
that is, in essence, through the source directivity diagram. 
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4 Alternative approach in solving the REB equa- 
tion 



The REB equation solving method stated in the previous section is based on rep- 
resentation of quantity P(z,p,Lu) as a sum of the singular part P {z,p,uj) and 
the regular one X{zq,p,uj) . Regularity of the X{zq, p,uj) has allowed one to use 
the expansion into the Taylor series over variables p and uj at the equation (|lO|) 
right side and to transform the integral-differential equation ( p^ ) into the first order 
partial derivative differential equation (pT|). 

However, the stated approach is not the only possible. The REB equation can 
be transformed right away using Fourier-representation of the function P(z, p, cj) : 

P(z,p,w) = -^^^^2 y"d^qP(z:q:'^)exp(iqp) (16) 
Substitution of ( p6| ) into (^) gives the following equation for quantity P{z,p,uj) : 

^P{z,q,uj) ^ J dLj'Q{z;Lj;Lj') ^P{z,q,u;')exp (j,q${z;Lj;Lj')^ -P(z,g,w)| 

(17) 

The quantity P{z,q,uj) is a differentiable function even when P(z,p,uj) has 
some peculiarities. Therefore, in the invariant coordinate small-angle scattering 
approximation it is possible to use the following expansion: 

P(z, q, Lo') - P(z, g, u) + ^£(£lM(^' _ ^) . (18) 
Substituting (^8|) in ( p^ ) we obtain the partial derivative differential equation 



d ^ d 

—P{z,q,Lj)~ A {z,q,Lj)—P{z,q,Lj) - P{z,q,uj) S {z,q,Lj) = 0, (19) 



where 



S(z,q,w)= duj'Qiz' ;lu;lu') exp (^iq$(z'; w; w'), w)) 
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A (z, q,uj) ^ J da;'Q(z'; co; to') exp (^iq$(z'; lu; lu'), {uj' - uj) . 
The characteristic system 

dUJ ~ / ^ ^ dP ~ , ^ ^ ^ ^ , ^ 

— = - Aiz,q,uj), —=S{z,q,uj)P{z,q,uj) (20) 
dz dz 

with initial conditions P = Po((f, w), a; = wq at z = has the following solution: 

V{z,p,uj) ^ j^-^ J d^qPo(q,Wo)exp|iqp-|-^ dz' S [z', q, cj(z', q, wq] | (21) 

This solution of the REB equation turns into the expression ( |l^ for P when 

A {z,q,uj) — > . But the latter limit transition corresponds to the invariant 
coordinate small-angle scattering approximation used in the previous section under 
derivation of ( [T^ ) and subsequent expressions. Let us note, however, that in (^l|), 
in contrast with (^), any additional terms do not appear. It allows one to assume 
that in used approximation the ratio 
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X(z,p,c^) «P(z,p,c^) (22) 

is fulfilled. Additional arguments to the benefit of this assumption will be presented 
in the following section. 

5 Analysis of the solution of the REB equation 

We shall show, first of all, that the obtained solution satisfies to the energy conser- 
vation law. For this purpose it is necessary to carry out integration of the left and 
right sides of (|l5|) over lo and/? multiplied them previously by sin cos 0. The area 
of integration over angles is defined by the condition that both wave uj and wave 
Lo' achieve the same level z (since at level z their mutual scattering occurs). To 
satisfy this condition one should add factors 8 [hr{oj) — z] and Q [hr{uj') — z] to the 
integrand expression, where Q{x) is Heviside step function, hr{uj) is the maximum 
height which can be reached by a ray with parameters 9, ip . Now integration can 
be expanded over all possible values of angles, i.e., over interval ~ tt/2 for 6 and 
over interval 27r for (p . Then, (fst) becomes 



/ P(cj) sin6'cos6'(i(jj = J Pq{uj) sin 9 cos 6 du!+ 

dz' J duj J du'e [h,{uj) - z'] e [h,{uj') - z'] Q {z';uj,uj') [Po{uj') ~ Po{uj)] 

where P{uj), Po{ll!) is a result of integration of P{zq, p,u!) and Po(p,Lu) corre- 
spondingly over variable p . 

Due to antisymmetry of the integrand expression with relation to rearrangement 
of shaded and not shaded variables, the last term in (??) is equal to zero. Thus, 
equation (??) reduces to 

J P{zq, p,ui) sva 9 cos 9 dujd'^ p — j Po{p,uj) sin 9 cos 9 dajd^ p (23) 

expressing the energy conservation law: the radiation energy full flux through the 
plane remains constant regardless of scattering, as it should be in case of real 
(dissipative) absorption absence. It is not difficult to check that parity (^) is valid 
for the exact solution in the form (^) and also for the solution in the form (^). 

With relation to the solution in the form (9) the carried out discussion discovers 
one curious peculiarity. It appears that the radiation energy complete flux through 

the base plane is determined by the first term . The second one (X) gives 
zero contribution to the energy complete flux. 

Let us investigate in more detail the structure of quantity X{z, p, lo) in the in- 
variant coordinate small-angle scattering approximation. Proceeding to the Fourier- 
representation in the expression (^) produces 

X{zo,q,iu) = J^" dz'J duo'Q{z'-oj,u') \p {z\q,Lu')- P {z' .,q,oj) 



exp < iq 



^z-LO-uo')+D{zo,0,Uj)\] 



Employing regularity of function P (z, q, uj), the last expression can be written as 

/■^o d p (z,q,ui) ~ \ ^ 

X{zQ,q,uj) = j dz' A (2:',g,w)exp iqD{zo, z' ,uj) 

where quantity A {z,q,uj) is defined by ([l^). Thus, it becomes evident that limit 
transition A {z,q,uj) — > entails also X(zo,p,uj) — > . This property has 
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been established in section 4 with the aid of comparison of two variants of the 
REB equation solution. Now we see that its presence is determined by structure of 
quantity X(zo, p, w). 

Results of the present section give the weighty ground to believe that the radi- 
ation spatial-and-angular distribution is determined basically by the first term in 
the solution (^). The second term represents the amendment to the solution which 
can be neglected in the invariant coordinate small-angle scattering approximation. 
This statement validity can be checked under detailed research of properties of the 
obtained REB equation approximated solutions by numerical methods. 

6 Conclusion 

In the present work the heuristic basis for use of the invariant coordinate small-angle 
scattering approximation is considered under solving of the RTE for a magnetized 
plasma layer. Within the framework of this approximation two versions of the ana- 
lytical solution have been obtained. They describe spatial-and-angular distribution 
of radiation reflected from a monotonous plasma layer with small-scale irregularities. 

The final physical conclusions about influence of the multiple scattering effects in 
a layer of plasma on the spatial-and-angular characteristics of radiation are possible 
on the basis of detailed numerical research of the obtained solutions. Such research 
is a subject of other our works. 
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